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ABSTRACT 

The results of experiments on the formation of molecular hydrogen on low density and high density 
amorphous ice surfaces are analyzed using a rate equation model. The activation energy barriers for 
the relevant diffusion and desorption processes are obtained. The more porous morphology of the 
low density ice gives rise to a broader spectrum of energy barriers compared to the high density ice. 
Inserting these parameters into the rate equation model under steady state conditions, we evaluate 
the production rate of molecular hydrogen on ice-coated interstellar dust grains. 
Subject headings: dust — ISM; abundances — ISM; molecules — molecular processes 



1. INTRODUCTION 

The formation of molecular hydrogen in the interstellar 
medium (ISM) is a process of fundamental importance in 
astrophysics because H2 helps the initial cooling of clouds 
during gravitational collapse and enters, either in neutral 
or ionized form, most reaction schemes that make o ther 
molecules l|Du]ev Williams igs4 lWi11ia,ms Iflfl?^. It 
was recognized long ago that H2 cannot form in the 
gas phase ef Hciently enough to acco unt for its observed 
abundance ijGould & Salneter 196,1). It was proposed 
that molecular hydrogen formation takes place on dust 
grains that act as catalysts allowing the protomolecule 
to quickly release the 4.5 eV of excess energy in a time 
comparable to the vibration period of the highly vibra- 
tionally excited state in which it is formed. 

The process of H2 formation on grains can be broken 
up into a few steps as follows. An H atom approaching 
the surface of a grain has a probability ^ to become 
adsorbed. The adsorbed H atom (adatom) resides on the 
surface for an average time in (residence time) before 
it desorbs. In the Langmuir-Hinshelwood mechanism, 
the adsorbed H atoms quickly equilibrate with the 
surface and diffuse on the surface of the grain either by 
thermal activation or tunneling. Furthermore, atoms 
that are deposited on top of already adsorbed atoms are 
rejected. When two adsorbed H atoms encounter each 
other, a n H2 molecule ma y form with a certain prob- 
ability llWi11ia,ms Iflfi^ lHn11en ba,ch fc Salneter TflTf ; 
Hollenbach & Saln eter 197lt Hollenb ach et al. 1971 ; 
^moluchowski 1981; Eronowitz & Chane 198^; 

Dulev & Williams IQS d ipTr ronello & Averna 1988; 
Sjindford & Allamaiido lla, 199^ ^ ITakahashi et al. 1999 : 
Farebrother et a l. 200Q^. 

The sticking probability of H atoms on dust grains 
was calculated using semi-classical methods. Quantum 
mechanical calculations of the mobility of H atoms 
on a polycrystalline water ice surface, showed that 



tunneling between adsorption sites, would provide H 
adatoms with sufficient mobility to guarantee efficient 
H2 formation even at 10 K (Hollenbach & Salneter 19^ 
iHollenbach Salneter 1971t iHollenbach et al. 197f. 
The steady state production rate of molecular hydrogen, 
i?H2 (cm^^ s^^) was expressed by 

i?H2 = ^"-hWhctt^: (1) 

where nn (cm~^) and vyi (cm s~^) are the number den- 
sity and the speed of H atoms in the gas phase, respec- 
tively, (T (cm^) is the average cross-sectional area of a 
grain and ng (cm~^) is the number density of dust grains. 
The parameter 7 is the fraction of H atoms striking the 
grain that eventually form a molecule, namely 7 = ^77, 
where 77 is the probability that an H adatom on the sur- 
face will recombine with another H atom to form H2. 
The probability ^ for an H atom to become adsorbed on 
a grain surface covered by an ice mantle has been cal- 
culated by Buch and Zhang (1991) and Masuda et al. 
(1998). They found that ^ depends on the surface tem- 
perature and on the energy of the irradiation beam. For 
a surface at 10 K and beam temperature of 350 K, Ma- 
suda et al. (1998) obtained a sticking coefficient around 
0.5. 

Alternative mechanisms, which were shown to oper- 
ate in selected experiments on well characterized crys- 
talline surfaces, are the Eley-Rideal mechanism and the 
hot-atom mechanism. In the Eley-Rideal mechanism, 
an atom from the gas phase impinges on another atom, 
which is already adsorbed, and reacts with it before be- 
coming thermally accommodated with the surface. This 
mechanism has been verified experimentally for D atoms 
on a Cu surface (Re ttner & Auerbach 1996J . In the 
hot-atom mechanism ijHarris fc Ka.semo 19811) . an atom 
lands on the surface and moves from site to site at super- 
thermal speed until it reacts with another atom (before 
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it thermally accommodates with the surface). In this 
case, the atom might fail to give up all of its incoming 
kinetic energy, or it might use some of the energy gained 
in becoming confined to the surface to move across it at 
super-thermal energy. The two mechanisms can be dis- 
tinguished experimentally by measuring the cross-section 
for the reaction o r the energy carrie d away from the just- 
formed molecule l|Zecho et al. 2002) . The kinetic energy, 
as well as the vibrational and rotational states of the des- 
orbed molecule, are expected to depend on the surface 
mechanism through which it was formed. 

In the last few years, we have conducted a series of 
experiments to study the formation of molecular hydro- 
gen on dust grain analogues und er conditions relevant to 
astrophysi cal environ ments ( Pirronello et al. 1997a^ 

[Pirronello et al. 19974 PirronellQ et al. 19m 

iRoser et al. 20021 iRoser et al. 20m The experi- 
ments used well collimated beams of hydrogen and 
deuterium atoms. The production of HD that occurs 
on the surface of a dust grain analogue was measured 
both during the irradiation with the beams and during a 
subsequent temperature programmed desorption (TPD) 
experiment. In this case, the temperature of the sample 
is raised quickly to either desorb particles (atoms and 
molecules) that got trapped on the surface or to enhance 
their diffusion and subsequent reaction and/or desorp- 
tion. In order to disentangle the process of diffusion 
from the one of desorption, separate experiments were 
carried out in which molecular species were irradiated 
on the sample and then were induced to desorb. 

In an earlier set of experiments, the forma- 
tion of molecular hydrogen was studied on sam- 
ples of poly crystalline olivi ne and amorphous car- 
bon rPirron ello et al. 1997at iPirronello et al. 1997bl: 
iPirronello et al. 199911 and the res ults w ere ana- 
lyzed using rate equa t ion models llKa,tz~| , al. 199^ 
ICazaux fc Tielens 2002t ICazanx Tielens 2004h. In 
this analysis the parameters of the rate equations were 
fitted to the experimental TPD curves. These parame- 
ters are the activation energy barriers for atomic hydro- 
gen diffusion and desorption, the barrier for molecular 
hydrogen desorption and the fraction of molecules that 
desorb upon recombination. Using the values of the pa- 
rameters that fit best the experimental results, the effi- 
ciency of hydrogen recombination on the polycrystalline 
olivine and amorphous carbon surfaces was calculated for 
interstellar conditions. By varying the temperature and 
fiux over the astrophysically relevant range, the domain 
in which there is non-negligible recombination efficiency 
was identified. It was found that the recombination ef- 
ficiency is highly temperature dependent. For each of 
the two samples there is a narrow window of high effi- 
ciency along the temperature axis, which shifts to higher 
temperatures as the fiux is increased. For the astrophys- 
ically relevant fiux range the high efficiency temperature 
range for polycrystalline olivine was found to be between 
7 — 9K, while for amorphous carbon it is between 12 — 16 
K. 

More recently, these studies were extended to different 
forms of water ice, namely high density ic e (HDI) and low 
density ice (LDI) (Manico et al. 2001; Ro ser et al. 20021: 
iHornekaer et al.200Hi) . In these experiments it was seen 
that the type of amorphous ice affects both the kinetics 
of molecular hydrogen formation and its efficiency. We 



then carried out a detailed analysis of the new data us- 
ing a rate equation model in order to obtain useful phys- 
ical parameters for the modeling of chemistry in dense 
interstellar clouds. Specifically, the parameters of the 
rate equation model are fitted to the experimental TPD 
curves. These parameters are the activation energy bar- 
riers for diffusion and desorption of hydrogen atoms and 
hydrogen molecules as well as the population ratios of 
different types of adsorption sites for molecular hydro- 
gen. Using the values of the parameters that fit best the 
experimental results, the efficiency of hydrogen recom- 
bination on the amorphous ice surfaces is obtained for 
values of the hydrogen flux and the surface temperature 
pertinent to a range of interstellar conditions. 

The paper is organized as follows. In Sections |2l and 
El we describe the experiments to be analyzed and their 
results. The rate equation models are introduced in Sec. 
31 Subsequent analysis and results are presented in Sec. 
El followed by a discussion in Sec. Eland a summary in 
Seed 

2. REVIEW OF EXPERIMENTAL METHODS 

The experimental apparatus and measurement tech- 
niques were most recently described in Vidali et al. 
(2004). Here we give a brief outline. The apparatus con- 
sists of an ultra-high vacuum (UHV) chamber pumped 
by a cryopump and a turbo-molecular pump (operating 
pressure in the low 10"^" Torr range). The sample is 
placed in the center of the UHV chamber and is mounted 
on a rotatable UHV compatible liquid heHum continuous 
flow cryostat. By varying the flow of liquid helium and 
by using a heater located behind the sample, it is possi- 
ble to change the temperature of the sample over a wide 
range. In the thermal desorption experiments, the tem- 
perature is raised quickly from around 10 K to about 
30 K. Higher temperatures are obtained at slower rates, 
as when the sample temperature is raised to ~ 90 K to 
change the structure of the ice (see below) or to des- 
orb the ice layers. The temperature is measured by an 
iron-gold/chromel thermocouple and a calibrated silicon 
diode placed in contact with the sample. 

Atoms are sent to the surface of the sample via two 
triple differentially pumped atomic beam lines aimed at 
this surface. Each has a radio-frequency cavity in which 
the molecules are dissociated, optionally cooled to ~ 200 
K by passing the atoms through a cold channel, and then 
injected into the line. Dissociation rates are typically in 
the 75 to 90% range, and are constant throughout a run. 
Estimated ffuxes are as low as 10^^ (atoms cm^^s"^) 
ijRoser et al. 2002^ . 

The experiment is done in two phases. First, beams 
of H and D atoms are sent onto the surface, at a con- 
stant irradiation rate Fq, for a given period of time to 
(from tens of seconds to tens of minutes) while the sur- 
face temperature is maintained at a constant value Tq. 
At this time any HD formed and released is detected by 
a quadrupole mass spectrometer mounted on a rotatable 
flange in the main UHV chamber. In the second phase, 
starting at to (the TPD phase), the irradiation is stopped 
and the sample temperature is quickly ramped at a rate 
b (~ 0.6 K/sec on average), and the HD signal is mea- 
sured. The time dependence of the surface temperature, 
T{t), during typical TPD experiments is shown in Fig. 
1, and is approximated by a piece wise-linear fit. In the 
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ideal case in which the heating rate is nearly constant, 
the time dependence of the flux and surface temperature 
can be approximated by 

F{t)^Fo; T{t)=To: < t < to{2) 

F{t)^0; T{t)=To + b{t-to) : t>t„. (3) 

The HDI sample is prepared by directly depositing a 
measured quantity of water vapor via a stainless steel 
capillary placed in the proximity of the sample (a pol- 
ished copper disk). Before deposition, the water under- 
goes repeated cycles of freezing and thawing to remove 
trapped gases. The sample is held at 10 K or lower dur- 
ing deposition and the deposition rate is 8 layers/sec, for 
a total thickness of the order of 1200 layers. This method 
of preparin g the ice sample should pro duce a high density 
ice, or HDI l|Jenniskens fc Blake 1994^ . The other phase, 
namely the low density amorphous ice or LDI, can be ob- 
tained from HDI by heating it. This transformation is 
gradual but irreversible and occurs over a broad tempera- 
ture range starting from about 38 K. In our experiments, 
LDI is prepared by gradually heating the HDI to 90 K 
and holding the temperature at that value for at least 5 
minutes. Then the sample is cooled down to 10 K. In 
the experiments analyzed here, the LDI sample was kept 
at a temperature between 9 and 11 K during the irradi- 
ation phase. During the TPD phase, care was exercised 
so that the sample didn't reach a temperature close to o r 
higher than 38 K l|Roser et al. 2nn2HVida]i et al. 20041). 
Low and high density amorphous ice were shown to have 
similar c haracteristics and structure as interstellar ic e 
mantles llMaver fc Pletzer 198(iH.Tenniskens et al. 199^. 
The amorphous ice structure, which contains many 
pores with diameters ranging between 15-20 A, can 
adsorb H2 molecules, that can diffuse into ice micro- 
pores where they c an be trapped ijMaver fc Pletzer 198(1 
ILangel et al. 199411 . 

3. EXPERIMENTAL RESULTS 

Irradiations with beams of H and D ("H+D" there- 
after) were done both on LDI and HDI in order to ex- 
plore the formation processes of HD molecules on these 
amorphous ice surfaces. In a separate set of experiments, 
beams of HD and D2 molecules were irradiated on the ice 
surfaces. The latter experiments do not involve forma- 
tion of molecules. However, they allow us to isolate and 
better analyze some of the parameters relevant to the 
molecular formation processes explored directly in the 
former experiments. They also support our hypothesis 
that molecules that consist of different hydrogen isotopes 
behave quite similarly (although not exactly the same) 
as HD molecules and all of them could be treated by the 
same general rate equation model. 

Using the methods described in the previous Section, 
the following experiments were carried out. The H+D ir- 
radiation runs were performed with different irradiation 
times (2, 4, 8, 12 and 18 minutes on LDI; 6, 8 and 18 
minutes on HDI). The surface temperatures during ir- 
radiation were Tq ~ 9.5 K in the LDI experiments and 
To ~ 14.5 K in the HDI experiments. In the experi- 
ments with molecular beams the irradiation time was of 
4 minutes. The surface temperatures during irradiation 
of HD and D2 on LDI were Tq ~ 9 K and Tq ~ 10 K, 



respectively. In Table 1 we present a list of the exper- 
imental runs analyzed in this paper. During the TPD 
runs, the sample temperature is monitored as a function 
of time. The time dependence of the sample temperature 
during typical temperature ramps is shown in Fig. 1 for 
LDI (solid line) and HDI (dashed line). The temper- 
ature ramps deviate from linearity but they are highly 
reproducible. The symbols show the experimental mea- 
surements and the lines are piecewise-linear flts. For the 
LDI experiments a good fit is obtained with two linear 
segments, with an average heating rate of 6 ~ 0.6 K/sec 
over the whole ramp. For the HDI experiments a good fit 
is obtained with three linear segments, with an average 
heating rate of 5 ~ 0.13 K/sec. Although the temper- 
ature ramps deviate from linearity, most of the struc- 
ture of the TPD peaks takes place within a temperature 
range in which it is linear to a very good approximation. 
The analysis of the data is done using these temperature 
ramps, thus any possible effect of the deviation from lin- 
earity is taken into account. Such deviations are found 
to be minor. 

The desorption rates of molecules vs. surface temper- 
ature during the TPD runs are shown in Figs. 2-7. In 
Fig. 2 we present the desorption rate of HD molecules 
after irradiation by HD and H+D on LDI. In both cases 
the TPD curves are broad. In the case of HD irradi- 
ation, three peaks can be identified in the TPD curve. 
For H+D irradiation there are only two peaks. Interest- 
ingly, they coincide on the temperature axis with two of 
the three peaks obtained for HD irradiation. In Fig. 3 
we show the desorption rate of HD molecules after ir- 
radiation by H+D on LDI for several irradiation times 
between 2 minutes and 18 minutes. As in Fig. 2, all 
these TPD curves exhibit two peaks. The position of the 
high temperature peak is found to be independent of the 
irradiation time, indicating that this peak exhibits first 
order kinetics. The low temperature peak shifts to the 
right as the irradiation time decreases, thus, this peak 
exhibits second order kinetics. Fig. 4 shows the desorp- 
tion rate of D2 molecules after D+D and D2 irradiation 
on LDI. These TPD curves are also broad, and are quali- 
tatively similar to those of HD desorption, shown in Fig. 
2. 

In Fig. 5 we show the desorption rate of HD molecules 
after irradiation by HD and H+D on HDI. Here the TPD 
curves are narrower and include only one peak. In Fig. 6 
we show the desorption rate of HD molecules after irra- 
diation by H+D on HDI for three irradiation times. The 
peak shifts to the right as the irradiation time decreases, 
namely, it exhibits second order desorption. Fig. 7 shows 
the desorption rate of D2 molecules after D2 irradiation 
on HDI. A single peak is observed, qualitatively similar 
to the one shown in Fig. 5 for HD. 

In the experiments in which beams of H+D atoms were 
irradiated on LDI and HDI samples, most of the HD 
molecules detected were formed during the heat pulse. 
Only a small fracti on of them were fo rmed during the 
irradiation process ijB.oser et al. 2002ll . This indicates 
that at least under our experimental conditions, prompt- 
reaction mechanisms iDulev fc Williams 19m or fast 
tunneling (Hollenbach et al. 197lll do not play a major 
role in the formation of molecules. 
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4. THE RATE EQUATION MODELS 

In the TPD curves studied here most of the adsorbed 
hydrogen is released well before a temperature of 40 K is 
reached. We thus conclude that the hydrogen atoms on 
the surface are trapped in physisorption potential wells 
and are only weakly adsorbed. The mechanism for the 
formation of H2 (as well as HD and D2) is assumed to be 
the Langmuir-Hinshelwood (LH) scheme. In this scheme, 
the rate of formation of molecular hydrogen is diffusion 
limited. This assumption is justified by the Langmuir- 
Hke behavior observed in the desorption curves. 

The analysis of thermal programmed desorption exper- 
iments usually starts with the Polanyi-Wigner expression 
for the desorption rate R(t): 



R{t) ^ vN{tY eM~Ed/kBT), 



(4) 



where N is the coverage of reactants on the surface, (3 
is the order of desorption, v is the attempt frequency, 
Ed is the effective activation energy for the dominant re- 
combination and desorption process and T = T{t) is the 
sample temperature. In the TPD experiment, first order 
(/? — 1) desorption curves R(t) exhibit asymmetric peaks 
with a sharp drop-off on the high temperature side. The 
position of the peak is insensitive to the initial coverage, 
determined by the irradiation time. Second order des- 
orption curves {f3 = 2) exhibit symmetric peak shapes. 
These peaks shift toward lower temperatures as the initial 
coverage is increased (Chan et al. 1978). An important 
assumption is that all activation energy barriers are cov- 
erage independent. This assumption may not apply at 
high coverage. However, at the low coverages obtained in 
the experiments analyzed here (up to ~ 1% of a layer), 
it is a reasonable assumption. 

We introduce two models for describing the TPD 
curves: a simplified model which succeeds in describing 
the main characteristics of the TPD curves and a more 
complete model which describes them using a more gen- 
eral and accurate view. 

4.1. The Simple Model 

Consider an experiment in which a fiux of H atoms 
is irradiated on the surface. H atoms that stick to the 
surface hop as random walkers. The hopping atoms may 
either encounter each other and form H2 molecules, or 
desorb from the surface. As the sample temperature is 
raised, both the hopping and desorption rates quickly 
increase. In the models used here, there is no distinction 
between the H and D atoms, namely the same diffusion 
and desorption barriers are used for both isotopes. 

In the models we assume a given density of adsorption 
sites on the surface. Each site can adsorb either an H 
atom or an H2 molecule. In terms of the adsorption of H 
atoms, all the adsorption sites are assumed to be identi- 
cal, where the energy barrier for H diffusion is and 
the barrier for desorption is E^^. As for the adsorption 
of H2 molecules, we assume that the adsorption sites may 
differ from each other. In particular, we assume that the 
population of adsorption sites is divided into J types, 
where a fraction fij of the sites belong to type j, where 
j = 1, . . . , J, and J2j A^j = 1- The energy barrier for des- 
orption of H2 molecules from an adsorption site of type 

Tides 
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This model is motivated by the TPD curves obtained 
after irradiation by HD and D2 molecules on LDI, shown 
in Figs. 2 and 4, respectively. These TPD curves are 
broad and can be divided into three peaks. We interpret 
this feature as an indication that there are three types 
of adsorption sites for molecules, which differ from each 
other in the energy barriers for desorption. 

Let iVn [in monolayers (ML)] be the coverage of H 
atoms on the surface, namely the fraction of adsorption 
sites that are occupied by H atoms. Similarly, let A^H2(i) 
(also in ML) be the coverage of H2 molecules that are 
trapped in adsorption sites of type j, where j — 1, . . . ,J. 
Clearly, this coverage is limited by the number of sites 
of type j and therefore iVH2(j) < Mj. Since we assume 
that each site can host only one atom or one molecule, 
the coverage cannot exceed a monolayer, and thus Nh + 

E,^H2(j)<l- 

Our analysis shows that for LDI one needs three types 
of molecular adsorption sites, namely J = 3, while for 
HDI good fits are obtained with J = 1. For the case of 
LDI we thus obtain the following set of rate equations 



N„ =F 



iVH2(2)=/X2a7VH' 



(5) 
(6) 
(7) 
(8) 



The first term on the right hand side of Eq. 10 rep- 
resents the incoming flux in the Langmuir kinetics. In 
this scheme H atoms deposited on top of H atoms or H2 
molecules already on the surface are rejected. F repre- 
sents an effective flux (in units of ML /sec), namely it 
already includes the possibility of a temperature depen- 
dent sticking coefficient. In practice we find that for the 
conditions studied here the Langmuir rejection term is 
negligible and it is thus ignored in the simulations. The 
second term in Eq. 10 represents the desorption of H 
atoms from the surface. The desorption coefficient is 



Wh = i^eM~E^"ykBT) 



(9) 



where v is the attempt rate (standardly taken to be 10^^ 
s~^), E^^ is the activation energy barrier for desorption 
of an H atom and T is the temperature. The third term 
in Eq. Q accounts for the depletion of the H population 
on the surface due to diffusion-mediated recombination 
into H2 molecules, where 



a = veyi^{-Ef^ /ksT) 



(10) 



J is E^-Hj). 



is the hopping rate of H atoms on the surface and E^^ 
is the activation energy barrier for H diffusion. Here we 
assume that there is no barrier for recombination. If such 
a barrier is considered, it can be introduced as discussed 
in Pirronello et al. (1997b, 1999). 

Eqs. describe the population of molecules on 

the surface. The first term on the right hand side of each 
of these three equations represents the formation of H2 
molecules that become adsorbed in a site of type j = 1, 
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2 or 3. The second term in Eqs. (j^J"® describes the 
desorption of H2 molecules from sites of type j, where 



WH2U) = i^exp(- 



ij)/kBT), 



(11) 



is the H2 desorption coefficient and E^^{j) is the activa- 
tion energy barrier for H2 desorption from an adsorption 
site of type j. The H2 production rate R is given by: 



^ = E^H.(jm(j). 



(12) 



This model can be considered as a generalization of the 
Polanyi-Wigner model [see Eq. Q]. It gives rise to a 
wider range of applicability, compared to Eq. In 
particular, it describes both first order and second or- 
der desorption kinetics (or a cor abination) for diffe rent 
regimes of temperature and fiux llBiham et al. 199^ . 

Experiments that involve irradiation by molecules are 
useful since they enable an independent evaluation of the 
parameters of molecular adsorption on the surface. Con- 
sider the case where a fiux of H2 molecules is irradiated 
on the surface. In this case Eq. ij^l is no longer relevant. 
The H2 molecules on the surface are distributed between 
the different types of adsorption sites, in proportion to 
the population ratios ^j. Each one of Eqs. lO-® can 
now be treated as independent and solved analytically. 
The TPD curve consists of three first order peaks, at tem- 
peratures determined by the desorption barriers E^^{j), 
j = 1,2,3. The sizes of these peaks are determined by 
the population ratios /ij . Given a TPD curve, the energy 
barriers E^^{j) can be obtained from the temperatures 
TmaxU) of the corresponding peaks. Solving for the max- 
imal desorption rate by taking its derivative with respect 
to the temperature we obtain 



ij) 



V exp 



rpdes 
H2 



ksTmaxij) 



(13) 



In general, the parameters fij, j — 1, . . . , J can be ob- 
tained by fitting the TPD curve as a sum of J Gaussians. 
Each of the /ij 's is obtained as the ratio between the area 
below the corresponding Gaussian and the total area be- 
low the TPD curve. In the case of HDI such procedure 
is not needed because the TPD curve can be fitted with 
a single type of molecular adsorption site. 

4.2. The Complete Model 

The simple model presented above provides good fits 
to the results of the experiments in which H+D, HD and 
D2 were irradiated on HDI (Figs. 5-7). In this case it was 
assumed that there is only a single type of molecular ad- 
sorption site. In experiments on LDI, where several types 
of molecular adsorption sites are assumed, the fitting of 
the TPD curves can be improved by using a more com- 
plete model. In this model the H2 molecules are allowed 
to diffuse between the different types of adsorption sites. 
For the case of LDI (J=3), the complete model takes the 
form 



F 



Nh, (1) = f^iaNa^ - an, {1)Nh, (1) 



(14) 



- W^H.(l)^H.(l) + MlE"H2 0')A^H.(j)) (15) 

J = l 

Nu, (2) = ^i2aNu^ - an, i2)Nu, (2) 

3 

- WnA2)Nn,{2) + ^i2J2an,{J)NnAJ)) (16) 

Nh, (3) = maNa^ - a^, (3)iVH, (3) 

3 

- iyH.(3)A^H.(3) + Ai3 5I"H.(j)AfH.(j)). (17) 

The first equation is identical to the corresponding one 
in the simple model, namely Eq. J^J. The first and third 
terms on the right hand side of Eqs. lfT5|l - ltr7jl are the 
same as those that appear in Eqs. ©-lIHl of the simple 
model. The second and fourth terms describe the diffu- 
sion of H2 molecules between adsorption sites of different 
types, where an^j) = vexp{-E^^ /ksT) and E^f is 
the activation energy for hopping of an H2 molecule out 
of a site of type j. The second terms account for the 
hopping of molecules out of adsorption sites of type j, 
while the fourth terms account for their re-distribution 
between adsorption sites of the three different types. 

5. ANALYSIS OF EXPERIMENTAL RESULTS 

In this Section we use the rate equation models to fit 
the experimental results, and thus obtain the parameters 
that describe the diffusion and desorption of hydrogen 
atoms and molecules on ice surfaces. We use a combi- 
nation of analytical and numerical tools. In particular, 
the numerical integration of the rate equation models 
is done using the Runge-Kutta stepper. The result of 
the integration is a set of TPD curves that are a func- 
tion of the chosen set of parameters, some of which can 
be determined analytically from the experimental TPD 
curves. The data obtained from each TPD run includes 
the time dependence of the flux F{t) and temperature 
T{t). The temperature T{t) is measured directly via 
a thermocouple. The flux F(t) (ML /sec ) is estimated 
as described elsewhere ijVidali et al. 199^ . An approx- 
imate value for F{t), in the required units of ML /sec, 
can be obtained by integrating the TPD spectra, gener- 
ating the total yield of the various experiments. The flux 
is then obtained from the exponential flt indicated by 
Langmuir kinetics. It is important to stress that this is a 
lower bound value for the flux, and this value is reached 
only if there is no desorption of H atoms during the TPD 
runs. In the analysis we assume that each of the energy 
barriers E^^{j), j 1, 2, 3, represents a Gaussian distri- 
bution of energy barriers centered around it. Although 
the Gaussian distribution provides much better fits, the 
important features of the curves depend only on the pa- 
rameters mentioned above. The Gaussian distribution is 
inserted into the equations by effectively splitting each 
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of the Eqs. Ijl5|l - lll7p into about ten equations, each with 
a different energy barrier for H2 desorption according to 
the Gaussian distribution around its central value. 

5.1. Analysis of Experiments with Irradiation of 
Molecules 

The TPD curves obtained after irradiation of HD and 
D2 on HDI are narrow and exhibit a single peak (Figs. 5 
and 7). In contrast, the TPD curves obtained after ir- 
radiation of those molecules on LDI are broad and three 
distinct peaks can be identified in them (Figs. 2 and 
4). Since no surface reactions take place in these experi- 
ments, these peaks are of first order kinetics. Thus, one 
can use Eq. (tl^ll to obtain the activation energy E^^{j) 
for the desorption of H2 atoms that are trapped in an ad- 
sorption site of type j from the surface. This calculation 
is done in the context of the simple model. 

Using the complete model in the case of LDI, with 
j — 1,2,3, does not change the positions of the three 
peaks. The complete model allows diffusion of molecules 
between different types of adsorption sites. As a result, 
there is a fiow of molecules from the relatively shallow ad- 
sorption sites towards the deeper ones. Thus, the relative 
sizes of the three sub-populations of molecules adsorbed 
on LDI become temperature dependent. 

The fitting curves for the HDI experiments with HD 
and D2 are shown in Figs. 5-7. Some uncertainty ex- 
ists in the value of the energy barrier for HD desorption 
from HDI. The value ^^^^^'(l) = 68.7 meV was obtained 
from the experiments in which the irradiation times were 
relatively short (Fig. 5). The experiments in which the 
irradiation time was long give rise to the somewhat lower 
value of E^l^'il) = 65.5 meV (Fig 6). Additional uncer- 
tainty arises due to the choice of the choice of the attempt 
frequency v. A deviation by an order of magnitude in its 
value would modify the resulting energy barriers by 3-4 
meV. 

In the case of LDI, where three types of adsorption 
sites can be identified, the situation is more complicated. 
The energy barriers for desorption of H2 molecules from 
these three types of adsorption sites are obtained from 
Eq. Ijl3|l . The population ratios ^j, j = 1,2,3 are ob- 
tained from fitting the TPD curves using the rate equa- 
tions. These TPD curves are fitted very well using the 
complete rate equation model (Figs. 2 and 4). Within 
this paremetrization, a simplifying assumption used here 
is that the hopping rate q;h2 (j) out of an adsorption site 
of type j is linearly proportional to the desorption rate 
W}ir^{j) from the same site. The experimental data is 
best fitted when the ratio between these two barriers is 
around 1, namely aH2(j) = W%2(j). Assuming the same 
attempt frequency for the two processes this indicates 
that E^f = E^l''. This is an unusual result, since typi- 
cally the barriers for desorption are higher than the bar- 
riers for diffusion. However, it may indicate that in the 
porous LDI the dominant diffusion process is, in fact, 
desorption and re-adsorption on the surface within the 
pores. Note that the energy barriers for desorption of 
molecules are obtained directly using Eq. ifT^ll . thus the 
fitting process is needed only for the diffusion barrier. An 
additional parameter is the width of the narrow Gaussian 
distribution of the energy barriers for each of the three 
types of adsorption sites. The energy barriers for desorp- 
tion of molecules from the different types of adsorption 



sites and the population ratios of these sites are presented 
in Tables 2 and 3, respectively. 

5.2. Analysis of Experiments with Irradiation of Atoms 

In the experiments with irradiation of H and D atoms 
(H+D for brevity) or of D atoms only (D+D), the ad- 
sorbed atoms diffuse on the surface and recombine. The 
experimental TPD curves for the HD and D2 molecules 
can be fitted using either the simple or the complete rate 
equation models. The parameters that characterize the 
behavior of hydrogen molecules on the surface, namely 
the desorption barriers E^^{j) and the population ratios 
/ij , are available from the experiments with molecular ir- 
radiation. Thus, the only fitting parameters that remain 
are the energy barriers for diffusion and desorption of 
hydrogen atoms, namely E^^ and E^^. In Figs. 2 and 
3 one observes that the two peaks in the TPD curves ob- 
tained after irradiation with H+D on LDI coincide with 
the two higher peaks obtained after irradiation with HD 
molecules. Similarly, the peak obtained after irradiation 
of H+D on HDI coincides with the peak obtained af- 
ter irradiation of HD molecules (Fig. 5). The location 
of this peak is determined by the desorption barriers of 
HD molecules, no matter whether these molecules were 
deposited or formed on the surface. 

In Fig. 3 we present five TPD curves for H+D irra- 
diation on LDI, for different irradiation times. In each 
curve one can identify two peaks. The location of the 
high temperature peak is independent of the irradiation 
time. This is an indication of first order kinetics, namely 
that the molecules desorbed at this stage were formed 
at a lower temperature. They remained on the surface 
until the temperature became sufficiently high for them 
to desorb. The low temperature peak in each curve in 
Fig. 3 shifts slightly to the right as the irradiation time 
is reduced. This is a feature of second order kinetics, 
indicating that within this peak the molecules were des- 
orbed upon formation. Our interpretation is that these 
molecules were formed in shallow adsorption sites and 
thus easily desorbed at the temperature in which they 
recombined. These shallow adsorption sites can be iden- 
tified with those that correspond to the leftmost among 
the three peaks in the molecular irradiation experiment 
shown in Fig. 2. The fact that the same set of parame- 
ters can fit all the curves in Fig. 3, obtained for different 
irradiation times, as well as all the other experiments on 
LDI, provides very strong evidence in favor of our mod- 
els. Similarly, a single set of parameter provides very 
good fits to all the experiments on HDI. These parame- 
ters are presented in Tables 2 and 3. 

5.3. Implications to Interstellar Chemistry 

Using the parameters obtained from the experiments 
we now calculate the recombination efficiency of H2 
molecules on ice surfaces under interstellar conditions. 
The recombination efficiency is defined as the fraction of 
hydrogen atoms adsorbed on the surface which come out 
as molecules. In Fig. 8 we present the recombination 
efficiency of H2 molecules vs. surface temperature for 
LDI under fiux of 0.73-10-^ ML s"^ This fiux is within 
the typical range for both diffuse and dense interstellar 
clouds. This particular value corresponds to a density of 
H atoms in the gas phase of 10 (atoms cm~^), gas tem- 
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perature of 100 K and 5 x 10^'^ ad sorption sites per cm^ 
on the surface ijBiham et al. 200Tll . 

A window of high efficiency is found between 11-16K. 
At higher temperatures atoms desorb from the surface 
before they have sufficient time to encounter each other. 
At lower temperatures diffusion is suppressed while the 
Langmuir rejection leads to saturation of the surface with 
immobile H atoms and recombination is suppressed. At 
these low temperatures, the steady-state coverage is high 
and the rate equation model may not apply. As al- 
ready suggested in Pirronello et al. (1999), mechanisms 
such as Eley-Rideal or diffusion by tunneling, which are 
not taken into account in the model, may become sig- 
nificant in this regime. Nevertheless, we can speculate 
that if the Langmuir rejection remains significant even 
at higher coverage, the trend in recombination efficiency 
shown should remain qualitatively correct. Luckily, such 
low temperatures are rarely of astrophysical interest. At 
high temperature atoms desorb from the surface before 
they have sufficient time to recombine, thus inhibiting 
hydrogen recombination at high temperature. For both 
high and low density ice the calculated asymptotic val- 
ues of coverage are very low and well within the regimes 
of experimentation and subsequent numerical simulation. 
Consequently, the relevance of the model is justified. Our 
results indicate that recombination efficiency of hydrogen 
on both HDI and LDI is high in the temperature range 
of ~14-20 K which is relevant to the interstellar clouds. 
Therefore these ice surfaces seem to be good candidates 
for interstellar grain components on which hydrogen may 
recombine with high efficiency. 

6. DISCUSSION 

One of the important features demonstrated by these 
experiments is the relevance of the surface morphology 
to the recombination rate of H2 molecules. Surfaces of 
amorphous ice are difficult to model due to the paucity 
of morphological information on either actual ISM ice 
coated grains or their laboratory analogues. This dif- 
ficulty gives rise to some uncertainty about the role of 
quantum effects. In principle, quantum tunneling of 
H atoms on the surface should be considered. How- 
ever, these experiments, much like earlier experiments 
on a morphous carbon and p olycrystalline olivine sam- 
ples ijPirronello et al. 1997bli . indicate that the mobility 
of the hydrogen atoms is very low at low temperatures. 
Thus, the effect of tunneling appears to be small. Cazaux 
and Tielens showed that quantum effects might be im- 
portant at high coverage. However, at the low coverages 
of H atoms on interstellar grains, tunneling effects are 
expec ted to be very small as can be seen in Fig. 6b in 
Ref. l|Ca,za,ux fc Tielens 2flf)4D . obtained for low irradia- 
tion (exposure) times. 

The effect of the ice morphology is best seen in the 
difference between the types of adsorption sites and 
their depth. The more porous structure of the LDI 
is refiected in the broad distribution of the activa- 
tion energies for desorption from different adsorption 
sites which spans the range between 40 to 60 meV 
ij.Tenniskens et al. 199,'^^ . In comparison, the energy bar- 
riers for desorption from HDI exhibit a narrow distri- 
bution around 68.7 meV. The energy barriers for dif- 
fusion of H atoms on the ice surfaces are found to be 
higher than those obtained earlier for more smooth sur- 



faces such as the polycrystalline olivine analyzed by Katz 
et al. (1999). Thus, the mobility of H atoms on the highly 
amorphous ice is relatively low, in qualitative agreement 
with the findings of Smoluchowski (1983). The energy 
barriers for diffusion and desorption of H atoms and 
molecules on amorphous ice were calculated using a de- 
tailed model of the amorphou s ice structure and the in- 
teractions at the atomic scale ("Buch fc C zerminski 19911 
Such & Zhang 1991; Hixson et al. 1992). The energy 
barriers for H desorption were found to be distributed 
between 26 to 57 meV, while the barriers for diffusion 
were found to be between 13 to 47 meV. Our results for 
LDI are in good agreement with these calculations, while 
for HDI our barriers are somewhat higher (see Tables 2 
and 3). For H2 molecules they find a distribution of 
activation energy barriers for desorption in the range be- 
tween 400 and 1200 K, namely 33 meV - 133 meV. Our 
three barriers for the desorption of HD molecules are in 
the range of 46 - 61 meV, namely around the center of the 
calculated distribution. Thus, our model captures the ef- 
fects of the amorphous ice structure and interactions, in 
agreement with detailed calculations at the atomic scale. 

In the experiments we use a combination of the H and 
D isotopes. In most cases the production rate of HD is 
measured, complemented by some measurements of D2 
formation. Experiments measuring H2 formation are dif- 
ficult to carry out. This is due to the existence of back- 
ground H2 which is the most abundant residual gas in the 
well baked ultra-high vacuum chamber. Even the small 
percentage of undissociated molecules that are sent with 
the atomic beam can infiuence the results. The use of 
hydrogen isotopes enables us to perform measurements 
with much more accuracy, and obtain the important 
physical parameters. The extrapolation of our results to 
processes that involve only H and H2 is non-trivial. How- 
ever, the experiments in which the formation of HD and 
D2 were analyzed, enable us to examine the isotopic dif- 
ferences. The analysis shows that variations between the 
energy barriers involved in the HD and D2 production 
processes do exist. Nonetheless, the two processes are 
qualitatively similar and the only difference is in a small 
increment in the energy barriers of the D2-formation pro- 
cess in comparison to the HD. We thus conclude that the 
energy barriers we obtain provide a good approximation 
to those involved in the production of H2 on ice surfaces. 

Previous rate equation models, used in analy- 
sis of hydrogen recom bina tion on carbon and 
olivine jKatz et al. 199(t ICazanx fc Tielens 20021: 
iCazaux fc Tielens 20041) . assumed that a significant 
fraction (denoted hy 1 — ^) of the hydrogen molecules, 
desorb from the surface upon formation (note that /i in 
these three papers had a different meaning than in the 
present paper). The rest of the molecules (namely, a 
fraction of them) remain on the surface and desorb 
thermally at a later stage. The origin of the parameter 
/i remained unclear, but it was necessary in order to fit 
the data. 

The model used in the present paper provides a phys- 
ical motivation for the parameter fi, based on a differ- 
ent interpretation of the fact that part of the molecules 
quickly desorb. It is based on the broad distribution of 
the energy barriers for desorption of hydrogen molecules 
which exhibits three types of adsorption sites. According 
to this model the recombined molecules equihbrate with 
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the surface and reside in adsorption sites until thermal 
desorption takes place. Molecules formed in the most 
shallow adsorption sites desorb quickly because by the 
time they are formed, the surface temperature is already 
high enough to activate their thermal desorption. The 
fraction, fii of the shallow adsorption sites among all the 
molecular adsorption sites is obtained from the experi- 
ments in which molecules are irradiated. 

As mentioned above, in the models used here it is as- 
sumed that hydrogen molecules do not desorb immedi- 
ately upon formation. Instead, they stay trapped in the 
adsorption sites or diffuse between them until thermal 
desorption takes place. Consequently, one needs to con- 
sider mechanisms for the dissipation of the excess en- 
ergy acquired from the recombination process in order 
to prevent prompt desorption. Although classical molec- 
ular dynamics simulations for H2 formation on amor- 
phous ice performed by Takahashi et al. (1999) give a 
large average kinetic energy (~530 meV) to the desorb- 
ing molecules, both Roser et. al. (2003) and indepen- 
dently Hornekaer et al. (2003) showed experimentally 
that this kinetic energy is much smaller (~3 meV), in 
agreement with our assumption that the excess energy is 
dissipated elsewhere. There are various mechanisms for 
efficient heat transfer from the molecule to the surface. 
One possible mechanism may be due to the very irregu- 
lar structure of the sample surfaces. Both LDI and HDI 
surfaces are disordered. Even if molecules are promptly 
released upon formation, they do not necessarily go di- 
rectly into the vacuum but they might undergo a mul- 
tiple series of collisions in which part of their energy is 
released to the solid. These multiple hits might lead to 
a subsequent re-adsorption of the molecule on the ice 
surface. Such a mechanism was already identified for H 
and D atoms impinging on an amorphous ice particle 
iBuch & Zhang 199]|. 

The experiments done by our group as well as by 
Hornekaer et al. have explored several types of ice mor- 
phology. The experiments analyzed here were done on 
amorphous ice of more than 1000 ML, where we have 
used two different ices, LDI and HDI (both amorphous 
and porous). The ex periments by Hornekaer et. al. 
ijHornekaer et al.2flfl3|1 were done on ice that was pre- 
pared using the same procedure used by our group to 
produce the HDI, with 2000 ML thickness. Their atom 
beam fluxes were ~ 10^^ (atoms cm~^ sec~^), namely 
about one order of magnitude higher than the fluxes 
used in our experiments. The range of exposure times 
was comparable. Hornekaer et al. investigated the ki- 
netics of HD formation and measured the efficiency of 
recombination and the energetics of the molecules re- 
leased from the ice layer after formation. The efficiency 
values they obtained are close to those obtained by our 
group on amorphous ice. The energy distribution of 
molecules formed showed that at least in porous amor- 
phous ice, molecules are thermalized by collisions with 
the walls of pores (where they are formed) before they 
emerge into the gas phase. Hornekaer et al. performed 
TPD experiments in which they irradiated H and D 
atoms either simultaneously or sequentially after wait- 
ing a delay time interval before dosing the other isotope. 
On porous ice the results they obtained are consistent 
with a recombination occurring quickly after atom dos- 
ing due to a high mobility of the adsorbed atoms even 



at temperature as low as 8 K. This high mobility was 
attributed either to quantum mechanical diffusion or to 
the so called hot atom mechanism, where thermal ac- 
tivation is not likely to play a significant role at this 
temperature. This conclusion is sensible in light of the 
high coverages of H and D atoms irradiated in their ex- 
periments, which required the adsorbed atoms to diffuse 
only short distances before encountering each other. As 
already suggested by Pirronello et al. (2004a, b) and 
Vidali et al. (2005), the hot atom mechanism may 
be able to provide the required mobility. In this case 
H and D atoms retain a good fraction of their gas 
phase kinetic energy during the accommodation process 
Puch & Zhang 1991; Takahashi & Uehara 2003)- This 
enables them to travel on the ice surface and inside its 
pores for several tens of Angstroms exploring several 
adsorption sites and recombining upon encountering al- 
ready adsorbed atoms. However, in our experiments, be- 
cause of the low coverage, the number of sites explored 
by the hot atom is not large enough to encounter an ad- 
sor bed atom and react with it with significant probabil- 
ity I Vida li et al. 20dBI) . Thus, there is no contradiction 
between the interpretations that molecular hydrogen for- 
mation is dominated by the hot atom mechanism at high 
coverage and by thermally activated mobility of the ad- 
sorbed H atoms at low coverage. 

In addition, Hornekaer et al. studied non-porous ice 
with only 20 ML. In this paper we have not analyzed 
such non-porous samples and therefore cannot provide 
any comparison in this case. On the other hand, we 
have used LDI for which Hornekaer et al. (2003) do not 
provide analogous results that enable a comparison. Our 
results and analysis show that on LDI the recombination 
of H atoms is thermally activated, and it is not efficient 
at low temperatures. 

In their analysis, Cazaux and Tielens (2004) consid- 
ered a rate equation model with a single energy barrier 
for molecular desorption from carbon and olivine sur- 
faces. In their model chemisorption sites play a role 
in the recombination process of hydrogen molecules and 
each hydrogen isotope is treated differently. Although 
chemisorption sites might play a role in the recombina- 
tion process, these sites do not seem to play an important 
role, at least not at the conditions explored in our exper- 
iments (low temperatures and low coverage as expected 
for interstellar dense cloud environments). However, the 
analysis of Cazaux and Tielens may provide a clue on 
H2 formation in photon dominated regions. We have 
not treated each atomic isotope differently but instead 
we used different effective averaged values for the experi- 
ments involving H and D atoms and those involving only 
D atoms. Consequently obtained a model that requires 
fewer fitting parameters. 

7. SUMMARY 

Experimental results on the formation of molecular hy- 
drogen on amorphous ice under conditions relevant to in- 
terstellar clouds were analyzed using rate equation mod- 
els. By analyzing the results of TPD experiments and 
fitting them to rate equation models, the essential pa- 
rameters of the process of molecular hydrogen formation 
on ice surfaces were obtained. These parameters include 
the activation energy barriers for diffusion and desorp- 
tion of hydrogen atoms and molecules on the ice sur- 
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face. While we identify only one type of adsorption site 
for hydrogen atoms, three types of adsorption sites are 
found for molecules on LDI, with different activation en- 
ergies. The parameters that determine what fraction of 
the molecular adsorption sites belong to each type are 
also found. Our model enables a unified description of 
several first and second order processes that involve irra- 
diation by either hydrogen atoms or molecules, all within 
the framework of a single model. 

The rate equation model allows us to extrapolate the 
production rate of hydrogen molecules from laboratory 
conditions to astrophysical conditions. It thus provides 
a quantitative evaluation of the efficiency of various 
ice surfaces as catalysts in the production of hydrogen 
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Table 1 

List of the TPD experiments on low and high density amorphous ice analyzed in this paper. 



Ice 
Type 


Atoms/Molecules 
Type 


Irradiation Time 
(minutes) 


Irradiation 
Temperature (K) 


LDI 


H 1 D atoms 


2 


~9.5 


LDI 


H + D atoms 


4 


~9.5 


LDI 


H + D atoms 


8 


~9.5 


LDI 


H + D atoms 


12 


~9.5 


LDI 


H + D atoms 


18 


~9.5 


HDI 


H + D atoms 


4 


~14.5 


IIDI 


H + D atoms 


6 


~14.5 


HDI 


H + D atoms 


8 


~14.5 


HDI 


H 1 D atoms 


18 


~14.5 


LDI 


D atoms 


4 


~10 


LDI 


HD molecules 


4 


~9 


HDI 


HD molecules 


4 


~15 


LDI 


D2 molecules 


4 


~10 


HDI 


D2 molecules 


4 


~15 
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Table 2 

The energy barriers obtained by the fitting of the TPD curves for low and high density ice surfaces. 

IS THE BARRIER FOR ATOMIC DIFFUSION, ^^h'^'' IS THE BARRIER FOR ATOMIC DESORPTION AND E^"^ (j) , j = 1,2,3 ARE THE 

BARRIERS FOR MOLECULAR DESORPTION FROM SITES OF TYPE j. 



Material 


Molecule 
Type 


Ej}^ (mey) 


fmeV; 


£it;(l)fmeV; 


Ef{:^{2)fmey) 


£^n°f(''^)('meV; 


Low Density 
Ice 


IID 
D2 


44.5 
41.0 


52.3 
45.5 


46.5 
40.7 


52.8 
53.3 


61.2 
65.5 


Higli Density 
Ice 


HD 
D2 


55 


62 


68.7 
72.0 





Table 3 

The PARAMETERS OF THE POPULATION RATIO OBTAINED BY THE FITTING OF THE TPD CURVES FOR LOW AND HIGH 
DENSITY ICE SURFACES, /il , /i2 AND /!;-{ ARE THE POPULATION RATIOS OF THE DIFFERENT ADSORPTION SITES OF THE 
MOLECULES AND (T IS THE STANDARD DEVIATION PARAMETER FOR THE ENERGY BARRIERS'S GAUSSIAN DISTRIBUTION 
AROUND THE CHARACTERISTICS ENERGY BARRIERS FOR MOLECULES DESORPTION. 



Material 


Molecule 


Ml 


M2 




<T (meV) 




Type 










Low Density 


IID 
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0.26 
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Fig. 1 

Fig. 1. — The time dependence of the surface temperature during typical temperature ramps in TPD experiments on LDI (steeper line) 
and HDI (less steep line). The symbols show the experimental measurements and the lines are piecewise linear fits. The irradiation phase 
is not shown. 
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Fig. 2 

Fig. 2. — TPD curves of IID desorption after irradiation with HD molecules (o) and II+D atoms (+) on low density ice. The irradiation 
time is 4 minutes. The solid lines are fits obtained by the complete rate equations model. 
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Fig. 3 

Fig. 3. — TPD curves of HD desorption after irradiation with H+D atoms on LDI. The irradiation times are 2 (*), 4 (o), 6 (□) 12 (+) 
and 18 (x) minutes. The solid lines are fits obtained by the complete rate equation model. 
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Fig. 4 

Fig. 4. — TPD curves of D2 desorption after irradiation with D2 molecules (o) and D+D atoms (+) on LDI, fitted by the complete 
rate equation model (solid lines). The irradiation time is 4 minutes. 
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Fig. 5 

Fig. o. — TPD curves of IID desorption after irradiation witli IID molecules (o) and H+D atoms (+) on high density ice, fitted by the 
complete rate equations model (solid lines). Irradiation time is 4 minutes. 
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Fig. 6 

Fig. 6. — TPD curves of HD desorption after irradiation with H+D atoms on HDI, fitted by the complete rate equations model (solid 
lines). The irradiation times are 6 (A), 8 (□) and 18 (o). 
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Fig. 7 

Fig. 7. — TPD curves of D2 desorption after irradiation with D2 molecules on HDI (o), fitted by the complete rate equations model 
(solid line). The irradiation time is 4 minutes. 
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Fig. 8. — Recombination efficiency of molecular hydrogen at steady state on LDI as a function of the temperature T (K), using the 
parameters obtained from experimental measurements of HD desorption curves. The fluxes of hydrogen atoms is 0.73-10~® ML s~^. 
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Fig. 9. — Recombination efficiency of molecular hydrogen at steady state on LDI vs. the temperature T (K), using the parameters 
obtained from experimental measurements of D2 desorption curves. The flux is the same as in Fig. 8. 
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Fig. 10. — Recombination efficiency of molecular hydrogen at steady state on HDI vs. the temperature T (K), using the parameters 
obtained from experimental measurements of HD desorption curves. The flux is the same as in Fig. 8. 



